On Analytic Properties of the Photon Polarization Function in a Background 

Magnetic Field 



W. F. Kao, Guey-Lin Lin and Jie-Jun Tseng 
Institute of Physics, National Chiao Tung University, Hsinchu 300, Taiwan 

(February 1, 2008) 

We examine the analytic properties of the photon polarization function in a background magnetic 
field, using the technique of inverse Mellin transform. The photon polarization function is first 
expressed as a power series of the photon energy ui with ui < 2m e . Based upon this energy expansion 
and the branch cut of the photon polarization function in the complex u) plane, we compute the 
absorptive part of the polarization function with the inverse Mellin transform. Our results are valid 
for arbitrary photon energies and magnetic-field strengths. The applications of our approach are 
briefly discussed. 
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The behavior of a charged particle in a background magnetic field is rather well known. The energy of the charged 
particle is quantized according to the Landau levels. For a non-relativistic electron moving in a uniform magnetic 
field along the +z direction, its energy levels are given by 

where lu c = eB/m e (e > 0) is the cyclotron frequency of the electron and s z is the electron spin projection along the 
+z direction. For a relativistic electron (positron), the energy quantization becomes 

= ml + P 2 z + eB(2n + l + 2s z ) , (2) 

with the following correspondent wave function 

V^,p B ,p,, Sz (7,t) = exp(~iE n t + ip y y + ip z z)F± p ^^(x'), (3) 

where F^ Pz §*( x ') ^ s a four-component spinor with x 1 = x — p y /eB. The detailed form for F^ Pz s ^ can be found, for 



example, in 

While the wave function of a charged fermion in the background magnetic field is well understood, it is often rather 
non-trivial to compute a physical process occurring in a background magnetic field, particularly, if there are more 
than one charged fermions involved in the reaction. A famous example is the pair production process 7 — > e + e~, 
which is relevant to the gamma-ray attenuation in the neutron star Q , and was first studied by Toll J3| and Klepikov 
independently. The approach by Toll and Klepikov is based upon a direct squaring of 7 — ► e + e~ matrix elements 
using exact electron and positron wave functions in a background magnetic field. The most updated calculation using 
this approach can be found in ||. To understand this problem better, however, one should demonstrate that the 
pair production width can also be obtained from the absorptive part of photon polarization function in a background 
magnetic field. In this regard, Tsai and Erber Q obtained the absorptive part of the one-loop photon polarization 
function in the asymptotic limit to 3> 2m e and B <C B c = m 2 e /e. Their result was shown || to agree with that of Toll 
and Klepikov. However, in the aforementioned asymptotic limit, the threshold behavior of the pair-production width 
is completely absent. Later on Shabad obtained the absorptive part (and the dispersive part as well) of one- loop 
photon polarization function for a general photon energy and magnetic-field strength. In that work, the threshold 
behavior of the pair-production width was worked out explicitly. We remark that Refs. |j|,f7| employed Schwinger's 
proper-time representation for the electron and positron Green's functions Q inside the photon polarization functions. 
To our knowledge, Ref. |Q is the first work which shows that the proper-time representation for the photon polarization 
function gives equivalent pair-production width to that obtained from squaring the 7 — > e + e~ amplitude directly. 
Unfortunately, the manipulations of Ref. [Q are rather involved and the substantial details of them were given in some 
other unpublished preprints ||. It is not very clear how one can generalize the approach of Ref. [Q to other processes. 

In this work, we will provide an alternative derivation of the pair-production width (or equivalently the absorptive 
part of the photon polarization function) from the proper-time representation of the photon polarization function. 
We shall also outline the procedure of obtaining the dispersive part of the polarization function, which is relevant 
to the photon index of refraction. It will be clear that our approach is very straightforward and physically intuitive. 
Furthermore, it is applicable to many processes occurring in a constant background electromagnetic field. In fact, 
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there are growing activities on computing the two-point and .three-point current correlation functions in a constant 
background field, using the string-inspired world- line formalisma. The final results of these calculations are expressed in 
terms of multiple integrals similar to the proper-time representations. Our approach will be very useful for investigating 
the analytic properties of these integrals, hence the analytic behaviors of various current correlation functions. 

To illustrate our approach, let us begin with the proper-time representation of photon polarization function Tl^ in 
a constant background magnetic field 113] : 



n ^(«)--7^ / ds J i d He"*[(9 2 V-*?^o 



/ .;, / 

(4tt) 2 

- (tf\9\\nv - q\\n<i\\v)N\\ + {q\g^u - q±^q±v)N ± ] 

-e- ism i(l-v 2 )(q 2 g^- qtt q u )}, (4) 

where q is the photon four-momentum with = (w,0, 0,q z ) and q^_ = (0, q Xl q y , 0) for a magnetic field in the +z 
direction, 

9 1 — v 2 9 cos(zu) — cos(z) 9 . . 



with z — eBs, and 



Mi = 



cos(zf) — i>cot(z) sin(zw) 



N\\ = - cot(z) (l-v 



sin(z) 

wsin(zu)\ cos(zv) 



sin(z) / sin(^) 



cos(zw) v cot(z) sm(zv) cos(zv) — cos(z) , , 

N± = + \\ , ' +2 V / 3 - V ; . 6 

sin(z) sm(z) sm' 5 (z) 

The two independent eigenmodes of the above polarization tensor arc ej^ and which are respectively parallel 
and perpendicular to the plane spanned by the photon momentum q and the magnetic field B. They obey the 
eigenvalue equations e^(—q 2 g t _ iu + Il^ejj = q 2 + Ily =0 and ^ A _{—q 2 g^ v + U^e^ = q 2 + H± = respectively with 

H\\ t ± = ejf ± ILj„ejj j_. It turns out that II|| and IT_l are proportional to N\\ and N± respectively. We shall not discuss 
the contribution by No since it does not correspond to an independent eigenmode. The absorptive part of H|| _i_ gives 
rise to the photon absorption coefficient (pair-production width) via the relation = ImII|| t ±/u). 
To study the analytic properties of Eq. (0), we employ the sum rule developed by us |j~3]]: 



1 / d n 



n! {d(oj 2 ) nILlx 



J 2 =0 7T 



dy l± -y^-U l± y-^\. (7) 



where Mil and M± are threshold energies for pair productions with Mi — q z = 4m 2 and Mj_ — qj = m 2 (l 



_ II 

yl + 2B/B C ) 2 ; y\\,±_ = M 2 j_/oj 2 . One notes that the absorptive part of II||_L(a; 2 ) vanishes for the range < uj 2 < 

M? j_ 1 14 1 . Therefore one can effectively set the integration range of Eq. (Q) as from y\\ t ± = to y\\,x — °°- 
Now, it is easily seen that the derivatives of IIm j_ at the zero energy are proportional to the Mellin transform of 
/ty,j_ • y» y 2 = K\\ t j_ • uj/M^x- Once the l.h.s. of Eq. (0) is calculated, the absorption coefficients or ImHj| : ± can 
be determined by the inverse Mellin transform. 

In order to compute the derivatives of Him at the zero energy, we note that a rotation of the integration contour 
s — > —is is permissible for a photon energy below the pair-production threshold. This rotation turns the oscillating 
trigonometric functions in IIim into more manageable hyperbolic functions. In the limit that B <C B Cl one can 
perform an asymptotic expansion in Bj B c , which gives 



1 / d n \ _ 2a e m 2 e fB 2 sm 2 0\" T(3n - l)r 2 (2n) /6n + l,3n + l 

n! \d(u} 2 ) nIll± ) ^=o ~ n 35 2 m 2 ) T(n)T{in) V 4n + 1 



(8) 



1 For an overview on recent developments, see [11001 
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where 9 is the angle between directions of the photon momentum and the magnetic field, and the neglected terms 
are higher order in B / B c . Applying the inverse Mellin transform, we obtain absorption coefficients k\\,± in agreement 
with previous results by Tsai and Erber 0. Such results are however not satisfactory. As one can easily see, apart 
from the trivial mass factor m 2 e r 2n , the first term on the R.H.S. of Eq. (ph grows to infinity as n increases, no matter 
how small the ratio B/B c is. This implies that the disregarded higher-order terms are in fact non-negligible for a 
sufficiently large n. As a result, the large moments of k\\± are not accurately determined by the first term of the above 
asymptotic expansion. Hence the threshold behaviors of k\\.± are not seen after applying the inverse Mellin transform. 
To generate correct threshold behaviors for absorption coefficients, we shall not expand in the magnetic-field strength 
even if B <C B c . We will obtain photon absorption coefficients valid for arbitrary photon energies and magnetic- field 
strengths. 

It is simpler to analyze II|| _l in the special case that q z = 0, i.e., q • B = 0. Without losing the generality, we may 
further choose g M — (ui, q x ,0,0). It is easy to show that II|| = — (a e uj 2 / 4ir)Ii\\ and II ^ = — (a e q 2 /4TT)Yl± with 

[•CO P + l 

H\\{u> 2 , q 2 .) = / dz dv exp[— isifio] Nu , 

n ± (w 2 , q 2 c ) = dz dv exp[-is<f> }N±. (9) 



Clearly, for the general case that <7 M = (uj, q x , 0, q z ) with q 2 = q 2 cot 2 the photon polarization functions, denoted 
as ft|| _l(w 2 , ql, 9), are given by II|| (uj 2 , q 2 , 9) = -(a e Lo 2 sin 2 9/4-k)TI\\ (q 2 , g 2 ), with q 2 = to 2 - q 2 = u 2 - q 2 x cot 2 9; 

and II_i_(w 2 , q 2 , 9) = — (a e q 2 ./'iir)Il±(q 2 , q 2 ). Let us now return to the special case. Rotating the integration contour 
s — ► —is and performing the energy expansion in to, we arrive at 

2 

n^io 2 iq 2 x ) = (A(a a ,(3 a )+B(a b ,(3 b ) + B(a b ,(3' b ) + C(a c ,(3 c ) + C(a' c ,(3' c )) , (10) 



and 



n ± (c 2 , ql) = ^ (A(a a ,/3 a ) + C(a c , &) + C(a' c , #) + D(a d , &)) , (11) 



where the functions A, B, C and D are defined as 



/2yi-l roc i- ! 

^ ■ ■ "b" (rT~orX 

,/2\n — I rCO rl 



i (to' ) n - 1 r°° 

A(a,[3) = Y, K ^ PP '(l' 2 )( ™)" V^TTT / dz / dv(l- v 2 ) n ' 1 (exp[azv] + exp[-azi;]) , (12) 



(i 



B{a,f3) = -^< ipp ,( g ' 2 )(— )"i-^ r J dz z n - 1 exp[-z/3] J dv (1 - v 2 ) n - l+1 (exp[azv] + exp[-azv}) , (13) 

i (lo' 2 y i ~ 1 r°° r 1 

C(a,(3) = -^^ w ,( g ' 2 )(_)"i_A_ J dz z n-l exp{-zP} J dv v(l - v 2 ) n - 1 (exp[azv] - exp[-cm;]) , (14) 
D(a,p) = -]Ti^ pp ,( 9 '2)(_)»^_ J dz z n - 1 cxp[-z/3] J dv (1 - v 2 )^ 1 (exp[azv] + exp[-azv}) , (15) 
with £ = EZo Elo E™=o E^=o> -' 2 = ^V4m 2 , & B/B c , q' 2 ee , 2 /4m 2 , 



,A ( 2(-1)'+p+p (2 g ' 2 )T(/ + m + l) B ,a,\ K A { n, 
Hmpp'Kl ) (l-p)\p\(l-p')\p'\T(m + l) ' lm PP' yq > 2 lrnpp ' {q ' 



KO („' 2 \-K^ (n' 2 \( l + m+l \ K° (n' 2 \ 8(-l)' +P+P ' (2q' 2 ) l r(l + m + 3)(/ + 1) 

«i mpP '[q )-«i mpP 'lq l + 1 j, «i mpP <lq ) - {l + 1 _ p)[p[{l - x _ pl)¥[{l - 2)r(m - 1} • W 

The actual arguments in the functions A, B, C , and D are given by a a = a c = a' c = p' — p+ 1, a b = a' b = a c \ = p' —p, 
Pa = p+p'+2m+l+l/B', f3 b = p a -l, P' b = p a + l, p c = p a , P' c = p a + 2, and d = /3 +l. Note that the various indices 



2 For convenience, we shall suppress the subscripts of a and (3 except in some special cases. 
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in the summation arise as follows: the index n comes from the photon-energy expansion, I arises from the binomial 
expansion of ((1 — v 2 )uj' 2 + 2(cosh(zw) — cosh(z))g' 2 /zsinh(z))", p and p' arises from writing (cosh(zv) — cosh(z))' as 
a sum of exponential functions, and m is due to expansions such as sinh~ z = 2 l cxp[— Iz] X)m=o ^m" 1 ^ 1 cxp[— 2mz]. 
It may appear at a first glance that our expansion on sinh - ' z and other similar terms are not convergent at z = 0. 
For z £ (0, e), one should expand sinh~ ( z in powers of z rather than in powers of exp(— Iz). Fortunately, as e — > 0, 
the result obtained with a careful treatment of sinh~' z does reduce to the one with sinh - ' z expanded naively. 

The computations of A, B, C, and D are rather similar in nature. For illustrations, we will go through the details 
of computing A(a,/3) and its relevant quantities. First of all, the integration over v can be carried out to give 

A(a,f3) = T,^ K ^ PP '(l' 2 )( W r ^ n J l + i ) J o dz z n - 1 exp[-z/3] Fi(n - 1 + -; -z 2 ), (17) 

where qF\ is the generalized hypergeometric function. We then perform the z integration which gives 

s 4 , ,os 1 x„ (uj' 2 ) n ~ l T(n-l + l) ,n-l + l n-l + 2 , 3 a 2 , 

A(a, 0) = E^(^)^(^)" ^ +lr ; n _, + i ; , ;» + W )- (18) 

For later conveniences, it is desirable to disentangle the indices n and I. We do this by replacing n — I + 1 with n and 
replacing the summation ^2 with the summation = Y^Li Ya^o E™=o 2»' p=o- This then gives rise to 

=-^2.^«w'(9 )(*) (^7) j^ry 2^-,—^+-,-^). (19) 

It is straightforward to compute the absorption coefficient by the inverse Mellin transform. We essentially invert the 
sum rule given by Eq. (Q). However, we should remark that Eq. (^) is derived under the assumption q 2 = with the 
angle between q and B kept general. Using this assumption, one has q 2 = to 2 sin 2 9 with 6 being the angle between q 
and B. It is then clear that there is no essential difference for analyzing the branch cut of Ily j_ either in the complex 
<7 2 -plane or in the complex w 2 -plane. This is the reason we have chosen lo 2 as the variable to set up the sum rule in 

Eq. (|7]). However, in the current case, q 2 is completely general. Hence it is more appropriate to study the analytic 
structures of H\\ t ± in the complex g 2 -plane. We generalize Eq. (Q) into 



-nil 



m; 2 " 



du\\ t ± -Uu ± • (wk||,j.) , (20) 



where m 2 = 4m 2 , = m 2 (l + \f\ + 2B') 2 , and U\\ t ± = mjj j_/q 2 - Here we use the notations ITj| and K|| _j_ to 
emphasize that the above sum rule holds for a general angle between q and B. We also note that the differentiations 
with respect to q 2 are understood as acting on il||(g 2 ,g 2 ) and Il±(q 2 ,q 2 ) containing inside IT|| and Tl± respectively. 
Taking the limit q - B = in the above sum rule, we can calculate the absorption coefficients Ku± through the 
following inverse Mellin transform: 



K[U (u; 2 , q l,B) = ^- I dsF lu± ( S , q 2 x ,B)(Q) 2s , (21) 



c 



where F^ ± (n, g 2 , B) = Tl^ ± \ Q2 - with w = u' for || polarization and w = 2uj' / (1 + \/l + 2B') for _L polariza- 

tion. The integral transform can be easily performed with the formula 



2 2 . 



2tt« Jc r (s + 5) 2 2 2 y/TT J 1 _ x + x 2 z 2 



One can derive this formula using the integral representation of hypergeometric function 3-^1 (fl, b; c; z) — 
[T(c)/(T(b)T(c - b))} ftdt ^(l - t) c - b - 1 (l - tz)- a and the fact that ^- \ c dsx~ s u s - x = S(x - u). Therefore, 
we arrive at 

K \\,A(a,l3) = 7T-, 1 > (23) 

2uj ^V(i-^) + (^) 2 
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where ^2" = EtoEm=oEp' p =q- To understand the structure of k^a, we rewrite the denominator of the above 
equation as 
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with (0 + a)B' = 1 + 2l\B' , and (0 — a)B' = 1 + 2l^B' . By a simple kinematic analysis, one can show that the 
R.H.S. of Eq. ( p4| ) is just |pi )Z |/m e a/ = |p2,z|/w e Ci/ where pi jZ and £>2,z are the z-direction momenta of e~ and e + 
respectively; while l\ and I2 are the Landau levels occupied by e~ and e + . The pair-production threshold corresponds 
to pi,z = j»2,z = 0. This threshold behavior is seen explicitly from the step function in the expression for K\\ t A- It 
is worthwhile to point out that the variables a's and 0's, appearing first in Eqs. ( |l0| ) and ([ll]), do have physical 
meanings. Taking A(a, 0) as an example, different values of a and imply a different radius of convergence for the 
infinite series in u>' (see Eq. (fT^)). The absorptive part of A(a, 0) emerges once the scaled photon energy, a/, becomes 
greater than the radius of convergence. This is reflected by the inverse Mellin transform given by Eq. (p2|). 



It is desirable to further simplify Eq. (|23|). We first perform the summation over the index I in 



Irapp' ' 



00 / 



x l T(l + m + l) 



£ xVr(p + m+ l)L»(-x), (25) 



^ ^ (I - p)\(l ~ p')\ 

l=0p,p'=0 y ^' v ^' p,p'=0 

where x = —2q' 2 /B' = — gr/(2e.B), p — max{p,p'}, p = min{p,p'}, and l^j^ m is the Laguerre polynomial. To derive 
this equation, we have used the relations iiq (g- 0; x) = e x \Fi(0 — a; 0; —x) and L"(x) = iF\{—n;a + l;x) x C]\ +a . 
At this juncture, the summation ^ in Eq. ( |23|) has already turned into X^ P '=o Em=0' To P rocee dj we recall that 
the actual values for a and in this case are a a = p' — p + 1 and a = p + p' + 2m + 1 + 1/ B\ which correspond to the 
Landau levels h = p' + m + 1 and I2 = P + m. Therefore the above summation is equivalent to X)z°^=i X^=o X)L=o 
with I = min{Ji — 1, By collecting the relevant terms in «n A(a a ,p a )i we can further sum over the index m: 

^ m T(l- m +l)ni- m + l)T{m + l) = W^T)^ 1 ^ (26) 
where I = max{Zi — 1, I2}. With the above summations over I and m, we arrive at 

w fefe v/(i-^) + (f^) 2 

with! 

T M,ft) = (-l) 1+rA x rA e x ^(Aa + 1) m (_a,)Lj* (-a), (28) 
ZiZ 2 v ; v ; T(X A +r A + l) Xa Xa 

where Ta = I — I = |z"i — ^2 — 1|, and = I = (Zi + Z2 — |^i — ^2 — 1| — The type B contributions to kii arise from 

B(ab,0b) and B(a' b ,0' b ), We write down the results without repeating the details: 



00 00 fvB 

a e m e B z x x ^ *-h 



where 



KM = (-l) r ^ B e* Ti ^^ B l [ 1} LZ (-x)L% (-x) + (X B A B - 1), (30) 



3 The absorption coefficient is not written in a symmetrized form for saving the space. Nevertheless, the symmetrization with 
respect to h and h can be easily done as one wishes. 
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with tb = \h — h\, As = (l-y + h ~ \h — h\)/%- We note that the contribution by B(a' b ,f3' b ) is incorporated by the 
replacement As ^ A^ — lin the above equation. Finally, the type C contributions are given by 



a e m e B 



/2 oc oo 



where 

(32) 

with rc = and Ac = A^. We like to emphasize that, in k\\, the relation between ui' 2 and q' 2 is still kept general. It 
is the dispersive part of Hii that determines the relation between uj' 2 and q' 2 . The dispersive part can be calculated 
using the Kramers- Kronig relation: 



- , ,9 ,9. P f°° Imri||(W 2 ) , 
Ren || (c' 2 ,g' 2 ) = - / dt \ (33) 



where P stands for evaluating the principle part of the integral, while the dimensionless function Ily (a/ 2 , q' 2 ) is related 
to Il||(w 2 ,g 2 ) by II||(a; 2 ,g 2 ) = — (a e o; 2 /47r)ri|| (uj' 2 , q' 2 ), as has been shown right above Eq. (||). It is well known that 



RelT is relevant to the photon refractive index n\\(u>') through the equation q 2 + RelF = 0. Hence 



n,| OS) = 1 + r dt . (34) 

For illustrations, we evaluate n\\ for a below threshold energy u> < 1 and a super-critical magnetic field B 3> B c - 
To the leading order in B c /B, it suffices to include the lowest Landau-level contribution on the R.H.S. of the above 
equation. Such a contribution is contained in K\\Bi i-e., 

h=i 2 =o a e m e B' 2q' 2 

= cxp , (35 

H' S 2a/ 2 yV 2 - 1 V £' V ' 

Combining Eqs. (p4[) and (pq), we find 



, ot e B' ( 1 „ i 
U||(w ) = 1 + -r — arctan \ — — 1 , (36) 

" V ' AllUJ ^^'2(1-^2) VI-CJ' 2 ' 

which agrees with the previous result Q. 

We now turn our attentions to the absorption coefficient k±. It is useful to rewrite Eq. ( pl| ) as 

K± (u;,qlB) = ^- f dsF ± (s,ql,B)(w') 2s , (37) 

where F±(n, q 2 , B) = ~j d7^p? ^-L L /2 =o- With written in this form, we can easily compute k± using the results 

from Comparing Eqs. ( |l0| ) and (pi]), we can easily show that kj_,a = — (q' 2 /uj' 2 )K\\ a, and k_l,c = — il' 2 / uj/2 ) K \\,c- 
Finally, it is slightly involved to compute the type D contribution. We obtain 



KID 



aemeB > r q> 2 ^^ T^ 2 (*)9 ((1 - + (f^) 2 ) 



with 

T£ i2 (x) = 8(-l) 1 +^^- 1 e-i^^L^+ 1 1 (-x)^- 1 (-x), (39) 
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with rr> = vb and Xo = Xb- Similar to the case of || polarization, the dispersive part of II j_ can also be calculated 
using the Kramers- Kronig relation. Since the technique is identical, we will not dwell upon this issue again. 

We have compared our absorption coefficients with those obtained by squaring the 7 — > e + e~ amplitude directlj 
Our results reduce to that of Daugherty and Harding in the special limit q 2 = which they have assumed. This 
is similar to what Shabad has demonstrated in Ref. |9) as he compared his result with that of Klepikov B in the 
above-mentioned limit for q 2 . Our approach differs from that of Refs. j7]Jj| in that Shabad performs the calculation 
in the beyond-threshold energy where the algebraic manipulations are rather involved and precautions are required, 
whereas we take the advantage of inverse Mellin transform which permits us to calculate the polarization function 
near the zero energy with a convenient energy expansion. We wish to stress again that our approach is ph ysic ally 
intuitive. We have written x as a multiple series in ui (photon energy), a and j3. It is possible to show that [p.5fl , for 
fixed values of a and /3, the series in w begins to diverge at the threshold energy uj t h = m e (\/l + 2l\B' + y/l + 2l 2 B') 
where the Landau levels l\ and l 2 are given by l\ = ((3 + a — l/B')/2 and l 2 = {(3 — a — l/B')/2. The divergence 
of the energy series beginning at lo = lo^ implies the existence of absorptive part beyond this point. We are able to 
compute the absorptive part as well as the dispersive part at any energy with the help of inverse Mellin transform 
and Kramers- Kronig relation. 

In closing, we have developed an integral-transform technique to compute the absorptive part of photon polar- 
ization function in a background magnetic field, while the dispersive part can be obtained via the Kramers- Kronig 
relation. Although we have chosen a special case q • B = for convenience, the result for a general angle be- 
tween q and B is easy to infer. Using the relations II|| (w 2 , q 2 , 8) — — (a e uj 2 sin 2 8/Att)YI^ (q 2 , q 2 ), YI±(uj 2 ,q 2 ,8) — 

— (a e q 2 ,/4ir)IL±(q 2 , <j 2 ), and k\^± = Imllm/w, it is clear that k\\ (uo 2 , q 2 , 8) = — (a e u> sin 2 #/47r)ImH|| (q 2 , q 2 ), and 
k^(u 2 ,q 2 ,8) = — (a e q 2 /47rw)Imn_i_(g 2 , q 2 )- For 8 = n/2, i.e., q • B = 0, K\\ — * K\\ = — (a e w/47r)Imil|| (to 2 , q 2 ) and 
k_\_ — ► k±_ = — (a e q 2 /47rw)lmllj_(a; 2 , q 2 ). Hence, for example, by comparing the above expressions for K|| and km, one 
realizes that K|| can be inferred from K|| by first dividing the latter by u>, then replacing the variable uj 2 by q 2 , and 

finally multiplying the entire expression by us sin 2 8. The procedure of obtaining k± from n± is also straightforward. 
There are other generalizations to the current work. For example, one may analyze the photon polarization function 
in a general background electromagnetic field, or study other current-current correlation functions under the same 
external condition. To our knowledge, the vector-axial vector correlation function relevant to weak interaction pro- 
cesses has not been analyzed as detailed as the fashion presented in this work. We shall report the results of such 
analysis as well as some technical details omitted in this letter in a forthcoming publication 
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